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Recently, different methods succeeded in calculating the spin dynamics at higher or- 
ders in the post-Newtonian (PN) approximation. This is an essential step toward the 
determination of more accurate templates for gravitational waves, to be used in future 
gravitational wave astronomy. We focus on the extension of the ADM canonical formal- 
ism to spinning binary black holes. Using the global Poincare invariance of asymptotically 
flat spacetimes as the most important guiding consistency condition, this extension can 
be constructed order by order in the PN approximation. We were able to reach a high 
order both in the spin power and the PN counting. 
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1. ADM Hamiltonian and global Poincare invariance 

In this article we consider asymptotically flat spacetimes and make use of standard 
(3+l)-splitting of spacetime into a stack of 3-dim. spacelike hypersurfaces. The 
Einstein equations are split into constraints and evolution equations. The speed of 
light c and the gravitational constant G are put equal to one. 

The ADM Hamiltonian 1 H[zl, P ai , S a{t)U) , hj^ , tt^ t ] basically is the ADM en- 
ergy expressed in terms of variables which possess standard canonical Poisson brack- 
ets (called canonical variables in the following). Here the matter variables z l a , P a i, 
and S a (i)(j) are the canonical position, momentum, and spin tensor of the a-th ob- 
ject, respectively, hj^ is the transverse-traceless part of the induced metric of the 
spacelike hypersurfaces and 7r^ T its canonical conjugate momentum (related to 
the extrinsic curvature). In order to calculate the ADM Hamiltonian, one only has 
to solve the field constraints (in the ADM gauge), not the evolution equations. 

The global Poincare algebra is a consequence of the asymptotic flatness and is 
represented by Poisson brackets of the corresponding conserved quantities P These 
quantities are the ADM energy H, total linear momentum Pi, total angular momen- 
tum tensor J^ , and the boost vector K%. The boosts have an explicit dependence 
on the time t and can be decomposed as Ki = Gi — tP{ . 



2. Hamiltonians linear in the single spin variables 

We treat black holes not as vacuum solutions of the Einstein equations. Instead we 
represent them by Tulczyjew's singular stress-energy tensor^ (in the covariant 
spin supplementary condition), which contains the matter variables. The matter 
variables are evolved by the Mathisson-Papapetrou equations.^ Then it is straight- 
forward to calculate the ADM energy or the other conserved quantities depending 
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on the matter variables appearing in T^ v , at least up to some order in the post- 
Newtonian (PN) approximation (using certain regularization techniques). However, 
the variables appearing in T^ v are, in general, not standard canonical and one has 
to redefine them. One approach is to find such a redefinition by requiring that the 
conserved quantities fulfill the global Poincare algebra. Then it must holcP 

Pi = p « - / d " x «%!? h ki;i > ( la ) 

a 

(lb) 

jJ'iTTiTT\ 1 f .3-. / ifcTT. TT _jArTTtTT\ 
~ x "can n kl,i) ~ / " 1 l~can n kj ~ "can "Asi J J 

which implies that a major part of the Poincare algebra is fulfilled. Requiring that P, 
and Jfj take this form indeed fixes the transition to canonical variables in a unique 
way at the next-to-leading order (NLO)P^ This allowed us to compute the NLO 
spin-orbit HamiltonianP which has already been derived earlier^ and agrees with 
non-canonical results P^l A new result was the complete NLO S1S2 Hamiltonianp^ 
which was later confirmed by a different method.^ Further, Gi was calculated and 
the full Poincare algebra was checked pi^i 

It was shown that at higher PN orders one also has to redefine the canonical 
momentum of the gravitational fieldPBy considering Eqs. ([1]), all variable redefini- 
tions necessary for the next-to-next-to-leading order can be found (up to a canonical 
transformation) M Variable redefinitions leading to canonical variables valid at all 
orders were recently constructed using an action approach!^ The used action basi- 
cally results from a minimal coupling of the flat space onef^ similar as in previous 
approaches.^ Next, all constraints, supplementary, gauge, and coordinate condi- 
tions are eliminated from the action. The variable redefinitions in question then 
transform this action into the form known to produce Hamilton's equations. 



3. Hamiltonians non-linear in the single spin variables 

If one knows the singular source terms of the constraints in terms of canonical 
variables, a straightforwardly calculated ADM energy will automatically be the 
sought-for Hamiltonian. This approach was used for higher orders in the spin for 
binary black holes. Source terms in canonical variables sufficient for H S 2 Sipi , H S 3 pi , 
H S 3 p2l H S 2g 2P2 , H S 2 S 2, H SlS 3, and Hs 2 sf were obtained from the Kerr metric in 
ADM coordinates.^ All these Hamiltonians are linear in G. The relation between 
Kerr parameter and Jij is crucial to obtain the canonical spin. Further, a leading 
order boost was applied to the Kerr metric, which allowed for leading corrections 
in the momenta. 

Another approach is to construct H and Gi depending on canonical variables 
directly from an ansatz and using the global Poincare algebra to fix the coefficients, 
with Pi and Jij still given by Eqs. ([I}. Ansatzes for H S 2 p2 , H S 2 p 2, Hg3 pi , H S 3 p2l 
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H S 2 g 2pi , H S 2 Sip2 , H S 4 , and H S 4 (again all linear in G) can indeed be fixed up to 
canonical transformation by {Gi, H} = Pi only. — Here also a corresponding ansatz 
for the source terms was used. The test mass case was used as a check. 

The static (P a j-independent) source part of the Hamilton constraint at NLO Sf is 
needed to fix the remaining degrees of freedom from the canonical transformation, as 
well as to get the missing G 2 -part of the NLO S 2 Hamiltonian (the linear-in-G part is 
given by H S 2 p 2). A general ansatz for these source terms, however, contains only four 
unknown constants.^ Two of these constants could be fixed by matching to the Kerr 
metric in ADM coordinates. The remaining two constants do not contribute to the 
Hamiltonian due to some cancellations. Finally, the complete NLO S 2 Hamiltonian 
can be calculated.^ The spin evolution agrees with an earlier resultP^D an( j th e 
test mass case was checked. If one wants to use the approach of the last section to 
derive the NLO Hamiltonian, one needs quadrupole corrections to Tulczyjew's 
stress-energy tensorPSl This, however, has not been tried yet. 

It should be noted that all Hamiltonians for binary black holes at the formal 2PN 
order are now known (for the formal counting, see, e.g., Ref. 1161 and also Appendix 
A of Ref.EJ). 
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